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Abstract We derive some formulas that rule the behaviour of finite differences under 
composition of functions with vector values and arguments. 



1 Introduction 

The Faa di Bruno formula [F] gives an expression for the n-th derivative (n > 1) of 
the composition fog of two functions / and g in terms of derivatives of / and g. 
Functions of class C°° will be called smooth. To simplify assumptions we shall always 
assume that / and g are both smooth functions. This formula has many versions, 
depending on the type of derivatives considered. Assume first that / and g are real 
valued functions of one real argument. Then the formula takes the following form 

where the sum is taken over all solutions (fci, . . . , 6„) G N" of the equation 

bi + 2b2 + . . . + nbn = n . 

We note that N = {0, 1, . . .} starts with zero, so that bi > for every i = 1, . . . ,n. 

Assume next that / and g are functions between euclidean spaces of possibly differ- 
ent dimensions, the domain of / containing the range of g. Let ei, . . . , be vectors in 
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the euclidean domain of g. The hst e = (ei, . . . , Cfc) will be referred as a fc-multi- vector. 
Given a = (ai, . . . , ak) G and a point x in the domain of g, we denote by D^^g^^e) 
the derivative o ■ ■ ■ o D^^g{x) of order |a| = cti + . . . + a^. Given a G N'^, let 

[a] = {/9GN'= : P<a} = {0,/3\P',...,P"'} . 
Then, Faa di Bruno's formula takes the form 

where the sum is taken over all solutions (61, ... , bm) G of the equation 

biP' + 62/3' + . . . + 6^/9™ = a , 

and a! = ai! ■ ■ ■ ttfc! when a = (ai, . . . , ak). Assume now the multi-index a belongs 
to the discrete cube I'' = {0, l}'^. In this case we have a\ = 1, /?M = 1, and bi G {0, 1}, 
which implies 6j! = 1, for alH = 1, . . . , m. Then formula ([2]) reduces to 

where the sum is taken over all solutions {hi, ... , bm) G of the equation 

bi(3^ + f>2/3' + . . . + b^r = « • 
We can simplify this formula a little more. 

Definition 1. Given a multi-index a E , we call partition of a with size r to any 
subset with r elements {a^ . . . , a'"} C — {0} such that + ■ ■ ■ + a"' = a. 

Of course the size r must range between 1 and \a\, since the smallest partition {a} 
has size 1, and > 1 for each a* in some partition. Then 

^°(/ ° gUe) = Yl ^"/^w • • • , D^^gAe)) (3) 

a^-\ \-a'^=a 

where the sum is taken over all partitions of a with size ranging from 1 to |a|. For- 
mula ([3]) is at the same time a special case and an extension of ([2]). Given a G N^, 
take e to be the r-multi-vector with components equal to Cj, for each i = 1, . . . , k, 
where r = \a\. Set a E P' = {0, 1}'' to be the multi-index with all components equal 
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to 1. Then -D"(/ o g)x{e) = D"^{f o g)x{e). Grouping and counting equal terms we 
can derive ([2]) from ([3]). 

After having written this article we found the reference |FL] where a discrete version 
of this formula is derived in terms of divided differences. Here we deal with finite 
differences in stead of divided differences, which means our formula holds for general 
functions of several variables. The history Finite Difference Calculus goes back a long 
way, parallel to that of Infinitesimal Calculus. We refer to [J] for a modern treatment 
of this calculus. We briefely recall some basic definitions in order to state our discrete 
version of Faa di Bruno's formula. Given a vector m G X, in some euclidean space X, 
let Tu be the translation ojoerator defined by {Tuf){x) = f{x + u). This operator acts on 
every space of functions / G . The difference operator along vector u, A„ : 
y^, is defined by A„ = id — r„, which corresponds to set (Au/)(x) = /(x + u) — f{x). 
Notice these operators always commute, i.e., A^oA^ = A^o A^i for all vectors u,v E X, 
since r^, o r„ = t^+v = t„ o r^. 

Definition 2. We call finite difference operator of order k to any composition of k dif- 
ference operators along possibly repeated vectors. Given a multi-vector u = {ui, ■ • • , Uk) G 
X'^ , we denote by A^ : Y^ Y^ the composite operator A^ = A„j o A„2 o ■ ■ ■ o A^^,. 
More generally, given a G N^, A" = (A^J^^o. . .o(Au^)"'° denotes a difference operator 
of order \a\. 

Next we introduce an algebra of symbolic finite difference expressions. Let x,y, . . . 
be symbols representing points, Ui, ^2, . . . be symbols representing vectors, and f,g,... 
be symbols representing functions. We denote by 2) the set of all symbolic finite 
difference expressions, which we define as the smallest set of expressions such that: 

1. if X is a symbol representing a point then x G 2); 

2. if -u is a symbol representing a vector then -u G D; 

3. if t G 2) and / is a symbol representing a function then f{t) G 2); 

4. if t, s G 2) then t + s G 2); 

5. if a G N'^, s,ti, . . . ,tk G 2) and / is a symbol representing a function then 

We consider as equal all terms which formally can be proved to be equal using property 
transformation rules of finite differences. Of course, depending on the interpretation 
given to the point, vector and function symbols, many terms in 2) will be formal but 
meaningless expressions. We define recursively the order of a term ord : 2) — > N: 
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1. if X is a symbol representing a point then ord(a;) = 0; 

2. if -u is a symbol representing a vector then ord('u) = 1; 

3. if t G D and / is a symbol representing a function then ord(/(t)) = 0; 

4. if t, s G 2) then ord(t + s) = min{ord(t), ord(s)}; 

5. if a G N^ s,ti,...,tk G D and / is a symbol representing a function then 
ord ^A^j^_ j = ai ord(ti) + . . . + ord(4). 

Given a term t = A"^^ tk)^^^^ ^ which is meaningful for some interpretation of its 
symbols (all functions being smooth), if all vectors Ui in t are small of order e then 
t is small of order eOrd(t) ^pj^g following theorem is our main result. Let X, Y and Z 
stand for linear spaces. 

Theorem A Given maps f G Y^ and g G X^ , and a multi-vector u G Z'' , 

KUo9){x)= A«^V.o....,AsW(^(^)) + ' (4) 

a^H \-a"=a 

where the ellipsis stand for higher order terms. 

Next theorem gives an explicit formula for (jll). 

Theorem B Given a G l'^, there are recursively defined sets, Aq, -A^i, . . . , asso- 
ciated with each partition ^ = {a^, . . . , a^} C of a, such that 

1. the sets ^^q^ ^q-i^ • • • ? ^ 

Q,r are pairwise disjoint; 

2. [3 eA^C [a], for /? = 0, a^, . . . , a"; 

3. I7I > \l3\, for every 7 G yi| - {(3}; and 

for any given maps f G Y^ and g G X^ , and any multi-vector u G Z'' , 

a:(/o^)(x)= Y1 (5) 

?={al,...,a'-}Gyc 

where tt| = X]7eyi* ^ufi'l^); /^'^ ^^^/i /3 = 0, q;\ . . . , a*". 

See theorem [H Note that Uq = g{x) + ■ ■ ■ , with remainder Xl^gyiS-io} ^udi^)-: 
for each z, = IS.'^ g{x) + with remainder Yl-yeA^ -^{0.^} ^Zdi^) ■ By 3. both 
these remainders are terms of higher order. 
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2 The Infinitesimal Formula 



Let D C X be an open domain. The tangent space T{D), and the tangent map 
Tf : T{D) T{Y), are defined to be T{D) ^ DxX, respectively Tf{x,u) = 
{f{x),Duf{x)). The chain rule shows that construct T is a functor, which essentially 
means that T{f o g) = Tf o Tg. Inductively, we can define higher order tangent spaces 
and tangent maps, by setting T^{D) = T{T'''^{D)) and T'=(/) = T{T^-^{f)). Then 
the iterated correspondence becomes also a functor. The tangent map of order 
fc, T^{f) can be explicitly expressed in terms of higher order directional derivatives 
of /. As we shall see the pattern of these expressions rules the behaviour of higher 
order derivatives under composition. To get an explicit expression for T^{f) we need 
some special notation to denote elements in T''\D). We shall call k-cuboid oi X to any 
family u = {uajaei'' of vectors in X indexed over the discrete cube I'' = {0, l}'^. Notice 
that any A;-cuboid u can be thought of as a pair of (A; — l)-cuboids, corresponding to 
restrict its indices to the two opposite faces {ak — 0} and {ak — 1} of the discrete 
cube l'^. Therefore, we can and shall identify the tangent space T''{X) with the set of 
all A;-cuboids of X: 



The /c-tangent space T^{D), to an open domain D, is the set of all A;-cuboids u e T^{X) 
whose base point uq belongs to D. The /c-tangent space over a point x E D \s the set 
T^{D) of all u e T^{D) such that uq — x. We shall use the multi-index derivative 
notation 



Because multi-indices are cumbersome to write, we shall adopt the following writing 
convention. Given a cuboid u e T^{X)^ Ui^^,„^i„ stands for the component where a 
is the multi-index (ai, . . . , ak) defined by aj = 1, if j G {zi, . . . , in}, «j = otherwise. 
Given a multi-index u = {ui, . . . , Uk) G X'' we shall denote by {{x; u)) the /c-cuboid w 
such that wq = x, Wi = Ui, . . ., = Uk, and Wi^^,,,^i^ = for all n > 2. With this 
notation is very easy to check that 

Proposition 1. If f E is a class C'^ map, then for all x E X, u E X^ , 



(X) = { M = {Ua)cel^ 



: X for all a e l'' 



}■ 



^:/(^) = p«.ro---o(^n,r/(x). 



rpk 



(/)(((^;«))) 



iD:f{x))^a'^ . 



The following kind of notation 





a^H \-a"=a 
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will always denote a sum taken over all partitions of a (see definition [T]) with size 
ranging from 1 to |a|. 

Proposition 2. Given a class map f G and a k-cuboid u G T^{X), writing 
fu = (T°' f u)a^jk , we have for each a & 

TVu= D:^....u.J{x) . (6) 

Q.1-1 \-a"=a 

Proof. This proposition is proved by induction. □ 
Let us call order of a term , „ f(x) to the sum la^l + . . . + la"!. Then, it is 

a-^ '*"' t/ \ / II II ' 

clear that expression (jS]) is homogeneous: all summands have order 
A couple of examples 

r(i'i'i)(/)(tl) = D„,,3/(x) + 

From proposition [H we see that the tangent map pattern ([6]) rules the behaviour 
of higher order derivatives under composition. 

Proposition 3. Given class maps f G and g G , and a multi-vector u G Z'^, 



qI-I \-a"=a 



Proof. 



D:{fog){x) = T^{f o g){{x-u)) 

= T\f)T\g){{x-u)) 



^D-'g{x),...,D^"g{x)-f(^) 

vl J L^n- 



Q H \-a"-=a 



□ 
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3 The Discrete Formula 



In order to characterize finite difference operators consider as before tlic discrete cube 
= {0, 1}*^ as a set of multi-indices. Given a E l'^ we write a = (cti, ■ ■ ■ , a^), wliere 
each ctj represents a binary digit, ai — or ai — 1. The set /'^ is partially ordered by 
the relation 

a < (3 ^ Oii < for all i = 1, . . . , /c . 

We also write \a\ — ai + . . . + ak and a-u — aiUi + . . . + akUk- A simple computation 
shows that 

Proposition 4. Given f e Y-^ , for all x E X, and u e X'^, 

A^/(x) = J](-l)H-l/3|/(a; + /3-ii) 

A key property of the difference operator A is the following kind of additivity. 
Proposition 5. Given f e Y-^ , for all x E X, and u,v E X, 

K+vf{x) = AJXx) + AJ{x + u) . (7) 

We define now a discrete equivalent of the fc-tangent map to / : X — > y. This will 
be a mapping Tk{f) : T^{X) T^{Y). The construct Tk will again be a functor. For 
that purpose we introduce the difference operator A : T^[X) T^[X) 

Au = (A'*M)^e/^ > where A"TI = . 

/3<a 

Notice that A^u — Uq. The operator A is invertible. Its inverse is the sum operator 
A-i : T'=(X) ^ 

The correspondence f eY^ ^ f^: T^{X) T^{Y), 

is obviously a functor. Thus, defining Tk{f ) = A o o A~^, the correspondence 
/ Tkif) is conjugated to / ^ Therefore, behaves functorially too. 
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Proposition 6. Given f G , for all x E X , and u G , 

Proof. It is enough to notice tliat A~^((a;; u)) = {x + a ■ u)^^jk. □ 

We shall say that any component Ua, of a fc-cuboid u, has order |a|. Then, we 
define recursively the order of a finite difference term A"^ n„/(^) to be the sum of 
the orders of terms Ui, . . ., Each term Ui can either be some cuboid component, as 
in proposition below, or else another finite difference term, as in theorem A. In both 
cases the formulas for Tafu and A"(/ o g){x), respectively, have a main part which is 
a sum of order \a\ terms, plus a remainder consisting of terms with order > |q;|. 

Proposition 7. Given f G andu G T!^{X), writing T^fu = {Tafu)a£ik, we have 
for each a E 

TJu= Yl K^....,u.J{^) + ■■■ (8) 

qI-I |~a"=Q: 

where the dots stand for a sum of higher order terms. 

Proof. This proposition follows from theorem [1], using the property ((71) to expand 
differences. □ 



Given a G we shall denote by Jq, the set of all partitions ^ = {a^, . . . , a""} of 
a. The multi-index obtained from a adding digit 1 at the end will be denoted by a<Ol- 
Therefore a<Ol belongs to I^^^ and has order |a<Ol| = |ct| + 1- Next lemma relates 
with ^aOi- 

Lema 1. Given ^ = {a^, . . . , a"} G 7a, consider the partitions of a()l 

= {001, a^OO, . . . , a"00} with size n + 1, 

= {a^OO, . . . , a^Ol, • • • , a"<>0} with size n, for 1 < i < n. 

These partitions exhaust CPaoi without repetitions, 

3^aOi = {^i '■ <i <n, and ^ G } . 

Given two /c-cuboids u, v E T'^^X), we shall denote by It; = ftl,vj the unique 
{k + l)-cuboid such that u^aoo = and Waoi = Va- 
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Lema 2. For all u, v e T''{X), 

1. A-i|lI,i;] = [A-in, /\-\u + v)]. 

2. A|ll,t;] = I An, At;- Am]. 

3. = IT,(/)1I, T^{f){u + v)-TkU)ul 



Lema 3. 

n 

'^w]ui,...,U„f{^) + ^Ui, ...,Ui^-i_,Vi, {Ui + i+Vi + i), {Un+Vn)f + W + Ui) 



i=l 



Proof. Follows by an iterated application of ([7]). □ 

We call maximum order of a partition ^ = {a^, . . . , a"} G 7^ to the number 
maxord(^) = max{|a-'^| , . . . , 

Theorem 1. Given a G , there are recursively defined sets, Aq, ^^i, • • • , A^^r asso- 
ciated with each partition ^ = {a^, . . . , a"^} C I'' of a, such that 

1. the sets ^^q, ^^^k^? • • • ? 

Q,r are pairwise disjoint; 

2. (3eA\c. [a], for /? = 0, a\ . . . , a^- 

3. < 7 < a and \^\ < maxord{^), for every 7 G /Iq — {0}; 

4- a* < 7 < a and \'~f\ < maxord{^) , for every 7 G A^^^ — {a*}; and 
for any given maps f G and g G , and any multi-vector u E Z'', 

To.fu= Yl (9) 

e={ai,...,a'}eTc 

where tt| = ^..^gy^? if7, /or eac/i /3 = 0, a^, . . . , a*". 
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Proof. This proposition is proved by induction in k. It is obvious when a E with 
k < 2. Assume it holds when k < n. Any given multi-index a G is of the form 
a()0, or a()l. Since the first case follows by induction hypothesis, we now restrict our 
attention to the second case. Given w G T''~^^{X), write it as a pair w = lu,v\, of 
/c-cuboids u, V E T^{X). Using lemma [3] we deduce 

Tc.oi{f){w) = T„oi(/)[M,tJl = T^U){u + v) - TM)u 

y A"e , fi4 + 4) - A" « fi4) 

^={Qi,...,Q"}eJ'<, 



^={al,...,a"}6J'c, 
n 



To finish the proof we just need to establish a one-to-one correspondence between 
summands above and partitions in CPq-oi. Using the notation introduced in lemma [H the 
partition ^0 is associated with the first summand, while the partitions (1 < i < are 
associated each with one of the subsequent n summands. Making these identifications 
we arrive at the equation between the previous sum and the following one: 

which, by lemma [1], is equal to 

€={/31,...,/3'-}gJ'<.oi 

We assume that = Xl^gyiS ""75 similar relations hold for tliI^q, etc. Then, 
matching terms, we arrive at the following recursive equations on the sets J^i: 



4^0 
•^000 


— -^0 


oO 


nlo 


— J\q 


ol 


/I Co 

a' 00 


= yif 

a 


.00 
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and 




= 4.00 

= 4,00 u 4,01 
4oo U AiOl U 4-00 



< 



if j = i 
if j <i 
if j > i 



where Aoi = {f3oi : f3 E A} for i = 0,1. These recursive relations ensure that 
Ta()i{f){w) has the correct development ([9]). It is now easy to check inductively that 



Theorem A in the introduction is a corollary of theorem [H Its proof is completly 
analogous to that of proposition [3] and, therefore, will be omitted. 

An algorithm that produces explicit expressions for Tafu has been devised, which 
was also used to confirm the correction of the above recursive definitions. Its source 
code can be retrieved from http:/ /ptmat.ptmat.fc.ul.pt/~pduarte/Research/FiniteDif- 
ferences/index.html. The following formulas were computer generated by this package. 



T{l,l)fiu) = Au.^J (Uo + U2+ Ml) + 



T{l,l,l)f(u) = Aui^2-J (^0 + U3+U2 + U2,3 + Ui+ Mi,3 + ^1,2) + 
^li,2,u,+U2,s+u,,J ("0 +U2 + Ui) + 



these sets satisfy the conditions 1., 2., 3. and 4. 



□ 



A3 




^lif ° 9)ix) = Aa2^ {g{x) + A,^g {x) + A^^g (x)) + 
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A^(/ o g){x) = Aa3^^„^ (gix) + A^^g (x) + A,.^g (x) + A^^g (x) + 

K,v,9 (x) + Al^^^g (x) + Al^^^g (x)) + 
AL,3W,A?„„3fl(x)/ i9{x) + A^^g (x) + A^^g {x)) + 

^\,^^g(x),A.,9{x)+A^.,,.,9{x)f (gi^) + ^v,9 {x) + A,,g {x)) + 
AA„j3(a;),A„2<;(x),A„3g(x)/ (f l^^)) 
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